Irreversible growth of binary mixtures on small-world networks 



o 
o 

(N 

< 



o 

B 



I 

O 

o 



> 

ON 

O 

o 



I 

o 
o 



X 
S3 



Julian Candia^ 

^ T/ie Abdus Salam International Centre for Theoretical Physics, 
Strada Costiera 11, 34OI4 Trieste, Italy 

Binary mixtures growing on small-world networks under far-from-equilibrium conditions are stud- 
ied by means of extensive Monte Carlo simulations. For any positive value of the shortcut fraction of 
the network (p > 0), the system undergoes a continuous order-disorder phase transition, while it is 
noncritical in the regular lattice limit {p — 0). Using finite-size scaling relations, the phase diagram 
is obtained in the thermodynamic limit and the critical exponents are evaluated. The small-world 
networks are thus shown to trigger criticality, a phenomenon analogous to similar observations 
reported recently in the investigation of equilibrium systems. 

PACS numbers: 05.70.Ln, 05.50.+q, 64.60.Cn, 75.30.Kz 



I. INTRODUCTION 

Complex networks are known to play a key role in 
the description of the structure and evolution of differ- 
ent mesoscopic and macroscopic systems. Indeed, the 
interacting parts of many natural and artificial systems 
can be interpreted as collections of linked nodes form- 
ing complex networks, whose structure and topology can 
be characterized in terms of statistical quantities such as 
their degree and path- length distributions, their connec- 
tivity, etc. 

Empirical observations show that the mean distance 
between a pair of nodes within a connected real network 
is in general surprisingly short (typically of a few degrees 
and only logarithmically dependent on the system size) , a 
phenomenon known as small-world effect. Moreover, the 
neighborhood of each node is observed to be, on average, 
highly interconnected 0. 

The well-studied classical random graphs, which are 
networks built by linking nodes at random, display the 
small-world effect but have much lower connectivities 
than usually observed in real networks. In this context, 
the small-world networks were proposed few years ago 
0, 0] as a realization of complex networks having short 
mean path-lengths (and hence showing the small-world 
effect) as well as large connectivities. Starting from a reg- 
ular lattice, a small-world network is built by randomly 
adding or rewiring a fraction p of the initial number of 
links. Even a small fraction of added or rewired links 
provides the shortcuts needed to produce the small-world 
effect, thus displaying a global behavior close to that of 
a random graph, while preserving locally the ordered, 
highly connected structure of a regular lattice. Indeed, it 
has been shown that this small-world regime is reached 
for any given disorder probability p > 0, provided only 
that the system size TV is large enough (i.e. TV > Nc, 
where the critical system size is A^c oc 1 /p) 0, Q . 

As a further step, recent works investigated the behav- 
ior of many standard models of Statistical Mechanics de- 
fined on small-world networks (see e.g. [MSt'! 8, 9, lOj), 
as well as on other classes of complex networks p] . In par- 
ticular, this was done for several equilibrium, Ising-type 



spin models. Besides their inherent theoretical interest, 
spin models defined on small- world networks can describe 
structural properties of polymer chains 0, 0, and 
have also been found useful in the study of social phe- 
nomena. For instance, spin states may denote different 
opinions or preferences, where the coupling constant de- 
scribes the convincing power between interacting individ- 
uals, which is in competition with the "free will" given 
by the thermal noise Moreover, a magnetic field 

can be used to add a bias that could be interpreted as 
"prejudice" or "stubbornness" [l5l |. 

Generally speaking, it was found that the structure 
and topology of the underlying complex networks affect 
dramatically the critical behavior of the models defined 
on them. For instance, it was found that the Ising model 
defined on a ID small-world network presents a second- 
order phase transition at a finite critical temperature Tc 
for any value of the rewiring probability p > 0, |^. 
Considering directed links, even the nature of the phase 
transition was found to change, switching from second 
to first order The ferromagnetic transition for the 

Ising model on small- world networks has also been stud- 
ied numerically by rewiring 2D and 3D regular lattices 

Much less attention, however, has been devoted so 
far to the investigation of nonequilibrium transitions on 
complex networks. Some simple nonequilibrium mod- 
els closely related to percolation were initially studied 
[T9I I20I I , while more recently a model for social in- 
teraction was investigated [2]| , in which the competition 
between dominance and spatial coexistence of different 
states in the nonequilibrium dynamics of Potts-like mod- 
els was examined. Within the context of these recent 
developments, the aim of this work is to investigate the ir- 
reversible growth of binary mixtures on small- world net- 
works. 

The growth of a binary mixture (or, adopting an equiv- 
alent magnetic language, a two-state magnetic system of 
up and down spins) can be studied by means of the so- 
called magnetic Eden model (MEM) [23, l23ll. a natural 
generalization of the classical Eden model |2J| in which 
the particles have an additional degree of freedom, the 
spin. In regular lattices, the MEM's growth process leads 
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to an Eden-like self-affine growing interface and a fractal 
cluster structure in the bulk, and displays a rich vari- 
ety of nonequilibrium phenomena, such as thermal order- 
disorder continuous phase transitions, spontaneous mag- 
netization reversals, as well as morphological, wetting, 
and corner wetting transitions. While the Eden model 
has extensively been applied to many different problems 
(namely, soot formation [25[. colloids [25l |. percolation 
|26j | , growth of cell colonies j2'j , crystal growth [2^ , 
etc), the MEM describes the aggregation of particles with 
a magnetic moment and can provide useful insight into 
kinetic phenomena of great experimental and theoreti- 
cal interest, such as the growth of metallic multilayers 
[snj l and thin films interacting with a substrate ) fluid 
adsorption on wedges > filling of templates imprinted 
with nanometer/micrometer-sized features |32o.M], etc. 
Moreover, interpreting the spin of the particles in a more 
general way, it can represent different atomic species in 
a binary alloy, impurities or defects in a growing crys- 
tal, the states of bacteria cells like Salmonella [33, the 
knowledge level of students in a classroom [s^, etc. 

The MEM growing on a small-world network could 
be considered, for instance, as representing the opin- 
ion spreading within a social group. According to the 
growth rules of the MEM, which are given in the next 
Section, the opinion or decision of an individual would 
be affected by those of their acquaintances, but opinion 
changes (analogous to spin flips in an Ising model) would 
not occur. However, as mentioned above, other interpre- 
tations of the model could also be possible in contexts as 
different as materials science, sociology, and biology. 

This paper is organized as follows: in Section 2, de- 
tails on the model definition and the simulation method 
are given; Section 3 is devoted to the presentation and 
discussion of the results, while the conclusions are finally 
stated in Section 4. 



II. THE MODEL AND THE SIMULATION 
METHOD 

In this work, we consider the one-dimensional, nearest- 
neighbor, adding-type small-world network model 
Is^ . Starting with a ring of N sites and N bonds, a 
network realization is built by adding new links connect- 
ing pairs of randomly chosen sites. For each bond in the 
original lattice, a shortcut is added with probability p. 
During this process, multiple connections between any 
pair of sites are avoided, as well as connections of a site 
to itself. Since the original lattice bonds are not rewired, 
the resulting network remains always connected in a sin- 
gle component. On average, pN shortcuts are added and 
the mean coordination number is (z) = 2(1 +p). Note 
that p can also be regarded as the mean shortcut fraction 
relative to the number of fixed lattice bonds. 

Once the network is created, a randomly chosen up 
or down spin is deposited on a random site. Then, the 
growth takes place by adding, one by one, further spins 



to the immediate neighborhood (the perimeter) of the 
growing cluster, taking into account the corresponding 
interaction energies. By analogy to the Ising model, the 
energy £^ of a configuration of spins is given by 

iv) 

where Si ^ ±1 indicates the orientation of the spin for 
each occupied site (labeled by the subindex i), J > 
is the ferromagnetic coupling constant between nearest- 
neighbor (NN) spins, and (ij) indicates that the sum- 
mation is taken over all pairs of occupied NN sites. As 
with other spin systems defined on complex networks, the 
magnetic interaction between any pair of spins is only 
present when a network bond connects their sites. 

Setting the Boltzmann constant equal to unity (fc^ = 
1) and measuring the absolute temperature T in units 
of J, the probability for a new spin to be added to the 
(already grown) cluster is defined as proportional to the 
Boltzmann factor exp{—AE/T), where AE is the total 
energy change involved. At each step, all perimeter sites 
have to be considered and the probabilities of adding a 
new (either up or down) spin to each site must be eval- 
uated. Using the Monte Carlo simulation method, all 
growth probabilities are first computed and normalized, 
and then the growing site and the orientation of the new 
spin are both determined by means of a pseudo-random 
number. Although the configuration energy of a MEM 
cluster, given by Eq.(^|, resembles the Ising Hamiltonian, 
it should be noticed that the MEM is a nonequilibrium 
model in which new spins are continuously added, while 
older spins remain frozen and are not allowed to fiip. The 
growth process naturally stops after the deposition of N 
particles, when the network becomes completely filled. 

For any given set of defining parameters (i.e. the net- 
work size N , the shortcut-adding probability p and the 
temperature T), ensemble averages were calculated over 
10^ different (randomly generated) networks, and consid- 
ering typically 50 different (randomly chosen) seeds for 
each network configuration. Since all normalized growth 
probabilities have to be recalculated at each deposition 
step, the resulting update algorithm is rather slow. In- 
volving a considerable computational effort, this work 
presents extensive Monte Carlo simulations that cover 
the whole shortcut-adding probability range < p < 1 
for different network sizes up to = 10''. 

III. RESULTS AND DISCUSSION 

The natural order parameter of a magnetic system is 
the total magnetization per site, i.e. 

which, in the context of this work, is to be measured 
on the completely filled network. However, since MEM 



3 





FIG. 1: Thermal dependence of the order parameter, for 
a fixed network size (A'^ — 100) and different values of the 
shortcut-adding probability p, as indicated. 



FIG. 2: Susceptibility as a function of the temperature, for 
different values of the shortcut-adding probability p and the 
fixed network size A'' — 100. 



clusters are grown from randomly chosen seeds, the en- 
semble average of the total magnetization is {M) ~ 0. 
Here, we will instead consider the absolute value of the 
total magnetization, \M\, as the order parameter, which 
is also appropriate to avoid spurious effects ari sing from 
finite-size spontaneous magnetization reversals |39l l40l | . 

Figure 1 shows plots of (|M|) as a function of T, for 
different values of p and a fixed network size, N = 100. 
The effect of increasing p at a fixed temperature is that 
of increasing the net magnetization (and, hence, the or- 
der) of the system. Indeed, larger shortcut fractions favor 
long-range ordering connections between distant clusters 
across the network. Considering instead a fixed value of 
p, we see that, at low temperatures, the system grows or- 
dered and the (absolute) magnetization is close to unity, 
while at higher temperatures the disorder sets on and 
the magnetization becomes reduced significantly. How- 
ever, fiuctuations due to the finite network size prevent 
the magnetization from becoming strictly zero above the 
critical temperature, and the transition between the low- 
temperature ordered phase and the high-temperature dis- 
ordered one becomes smoothed and rounded. 

Strictly speaking. Figure 1 is just showing evidence 
of pseudo-phase transitions, which might be precursors 
of true phase transitions taking place in the {N — > oo) 
thermodynamic limit. In the following, we will proceed to 
characterize in more detail this pseudo-critical state, by 
measuring other observables on finite-size systems. Fur- 
ther on, we will use standard finite-size scaling proce- 
dures to establish the phase diagram Tc vs p correspond- 
ing to the true phase transition in the thermodynamic 
limit, as well as to calculate critical exponents that de- 
scribe the behavior of the system at criticality. 

Let us now consider the magnetic susceptibility X; 
given by 

X^^{{M')~{\M\r). (3) 



For equilibrium systems, the susceptibility is related 
to order parameter fluctuations by the fluctuation- 
dissipation theorem. Although the validity of a 
fluctuation-dissipation relation in the case of a nonequi- 
librium system is less evident, we will assume Eq.Q to 
hold also for the MEM. Indeed, this definition of x proves 
very useful for exploring the critical behavior of this sys- 
tem, as shown in earlier studies of the MEM in regular 
lattices 1231 , as well as in other nonequilibrium spin mod- 
els mHa. 

Figure 2 shows plots of x vs T, for different values of 
p and a fixed network size, N = 100. As with the ther- 
mal dependence of the order parameter shown in Figure 
1, the order-disorder transitions signaled by the peaks 
of the susceptibility become rounded and shifted. For a 
given probability p and system size iV, we will define the 
effective finite-size "critical" temperature T^-^^{N;p) as 
the temperature corresponding to the peak of the sus- 
ceptibility. Although the transition temperature of finite 
systems is not uniquely and precisely defined, the suscep- 
tibility peaks become sharper as one considers larger sys- 
tems, and the effective finite-size "critical" temperatures 
tend to the true critical temperature in the (infinite-size) 
thermodynamic limit [23, |33, 113 • 

In the same vein, the heat capacity per site can be 
related to energy fluctuations as 

Cv = j^ {{E') - im . (4) 

Figure 3 shows plots of vs T corresponding to the 
same parameter values used before. Compared to the 
susceptibility, the heat capacity exhibits flatter shapes 
and broader peaks. 

Further insight can be gained by examining the nor- 
malized probability distribution of the magnetization, 
Pn{M). Fixing the network size {N = 1000), the be- 
havior of Pn{M) for different temperatures is shown in 
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FIG. 3: Heat capacity per site as a function of the tempera- 
ture, for a fixed network size (A'' = 100) and difi^erent values 
of the shortcut-adding probabihty p, as indicated. 



Figure 4 for (a) the regular lattice and (b) the small- 
world network with shortcut-adding probability p = 0.5. 
For high temperatures, in both cases the probability dis- 
tributions are Gaussian-shaped and centered at Af = 0, 
as expected for thermally disordered systems. However, 
their behavior is quite different at intermediate and low 
temperatures. 

For p = (see Figure 4(a)), the curvature of the dis- 
tribution is convex and exhibits a local maximum al- 
ways situated at M = 0, while other maxima develop 
at M = ±1. At very low temperatures, the curved shape 
flattens and the maxima at M = ±1 dominate the dis- 
tribution. Irrespective of the temperature, neither abso- 
lute nor local maxima arise at any intermediate values of 
the magnetization, i.e. different from M = (the fully 
disordered state) and M = ±1 (the completely ordered 
states). In fact, this noncritical behavior is in agree- 
ment with previous results for the MEM in the lattice: 
a cluster growing in a linear regular lattice from a sin- 
gle seed shows only a pseudo-phase transition with an 
effective "critical" temperature T^^^{N) that vanishes 
in the {N oo) thermodynamic limit '2^1 . An anal- 
ogous behavior was observed in the MEM grown in a 
stripped (l-t-lV-dimensional rectangular geometry using 
linear seeds |23|. 

In contrast, in the small- world network (see Figure 
4(b)) one observes the onset of two maxima located at 
M = ±Msp (0 < Msp < 1), which become sharper and 
approach M = ±1 as T is gradually decreased. The 
smooth shift of the distribution maxima across T ~ Tc, 
from AI = to the low-temperature nonzero sponta- 
neous magnetization M = ±Msp, is the signature of true 
thermally-driven continuous phase transitions l40| . 
Hence, critical behavior in the irreversible growth of 
MEM clusters arises from the presence of shortcuts in 
the small- world network. This result is a nonequilibrium 
realization of analogous phenomena observed in related 




FIG. 4: Normahzed probability distributions of the magneti- 
zation, for a fixed lattice size (A'^ = 1000) and different tem- 
peratures, as indicated, (a) The {p = 0) regular lattice, which 
exhibits local or absolute maxima at M = (the fully dis- 
ordered state) and M = ±1 (the completely ordered states). 
The sharp peaks at M = ±1 for T = 0.3 have been truncated, 
(b) The small-world network with shortcut-adding probabil- 
ity p — 0.5. The gradual onset of maxima at M = ±Msp 
(0 < Map < 1) across the transition, which become sharper 
and approach M = ±1 as T is decreased, is the hallmark of a 
true thermally-driven continuous phase transition. See more 
details in the text. 



equilibrium systems, as e.g. the Ising model in small- 
world networks generated from rewiring ID lattices 5, 
Notice, however, that this work focuses on nonequilib- 
rium processes that cannot be derived from the study 
of equilibrium systems. Hence, the phenomena observed 
here provide further and independent evidence on the 
effects of long-range interactions on dynamical systems 
which are defined on small- world networks. 

In order to explore further this phenomenon, we will 
extrapolate the finite-size "critical" temperatures to the 
thermodynamic limit and build the corresponding phase 
diagram, in which the critical temperature is given as a 
function of the shortcut-adding probability of the net- 
work. Moreover, this procedure will also allow to deter- 




FIG. 5: Effective transition temperatures T^''{N\p) for 
10^ < N < 10* and different values of p, as indicated (sym- 
bols). Fits to the data using the finite-size scaling relation, 
Eq.Q, are also shown (solid lines). See more details in the 
text. 



FIG. 6: Phase diagram T'c(p) vs p corresponding to the ther- 
modynamic limit, obtained using a finite-size scaling relation, 
Eq.Q. The MEM growing on small- world networks with any 
value of p > undergoes thermally-driven continuous phase 
transitions. In the regular lattice limit, p = 0, the system is 
noncritical. The dashed line is a guide to the eye. 



mine the critical exponents of the system. 

According to the finite-size scaling theory, developed 
for the treatment of finite-size effects at criticality and 
under equilibrium conditions [isl 0| , the difference be- 
tween the true p— dependent critical temperature, Tc(p), 
and the effective pseudo-critical one, T^^^ {N;p), is given 
by 

|^,(p)-T,^//(7V;p)|«7V-l/^ (5) 

where is the exponent that characterizes the divergence 
of the correlation length at criticality. 

The symbols in Figure 5 show the effective transition 
temperatures calculated for different network sizes in the 
range 10^ < N < 10'* and different shortcut fractions, 
as indicated. Recall that, as commented above, the size- 
dependent, pseudo-critical temperatures T^^^ {N]p) were 
determined for different system sizes and shortcut frac- 
tions from the maxima of the susceptibility (see Figure 
2). By means of least-squares fits of the finite-size scal- 
ing relation, Eq. (|SJ| , to these data one can determine both 
the true critical temperature of the system, Tdp), as well 
as the critical exponent v. The non-linear least-squares 
fitting procedure was implemented using the Levenberg- 
Marquardt minimization method |45l |. 

The solid lines in Figure 5 show least-squares fits to the 
data. Indeed, the finite-size scaling theory proves useful 
in describing the behavior of this nonequilibrium system 
near the critical region, since the scaling relation given 
by Eq.© provides an excellent fit to all the data. 

Figure 6 shows the phase diagram Tc{p) vs p, corre- 
sponding to the critical behavior of the system in the 
thermodynamic limit. The error bars reflect the statis- 
tical errors, which were determined from the fitting pro- 
cedure. As anticipated, for p > the system undergoes 
critical order-disorder phase transitions at finite critical 
temperatures: the small-world network geometry trig- 
gers criticality. Naturally, the global ordering imposed 
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FIG. 7: Plots showing the maxima of x for different network 
sizes in the range 10^ < N < 10* and different shortcut frac- 
tions (symbols). Also the corresponding finite-size scaling fits, 
given by Eq.lO, are shown for comparison (solid lines). 



by long-range shortcuts is weaker the lower the shortcut 
fraction, and hence Tc{p) decreases monotonically with 
p. The critical temperature vanishes for p = 0, which is 
the expected regular lattice limit behavior. 

As commented above, the same fits of Eq.© to the 
numerical data determine also the critical exponent i>. 
The obtained (p— independent) value is = 3.6 ± 0.4, 
where the error bar reflects the statistical error resulting 
from the fltting procedure. As in the case of equilibrium 
spin systems defined on small-world networks (see e.g. 
[3,^3), the universality class of this nonequilibrium sys- 
tem is not observed to depend on the shortcut density, 
provided that p > 0. 

An additional characterization of the critical behavior 



6 



of this system can be obtained by calculating the critical 
exponent 7, which describes the divergence of the suscep- 
tibility at the critical point. Using again the finite-size 
scaling theory 113,13 1 exponent ratio 'y/v can be re- 
lated to the peak of the susceptibility measured in finite 
samples of size N by 

Xma. CC N-^/r (6) 

The symbols in Figure 7 correspond to the maxima of x 
plotted against the network size for different values of p, 
as indicated, while the solid lines are fits to the data using 
this scaling relation, Eq.©. It turns out that ^/v — 
0.92 ± 0.04, where the error bar reflects the statistical 
error from the fit. Using this ratio and the value already 
obtained for we determine 7 = 3.3 ± 0.4. 

IV. CONCLUSIONS 

Binary mixtures growing irreversibly on small-world 
networks are studied numerically by means of exten- 
sive Monte Carlo simulations performed on the mag- 
netic Eden model. Firstly, evidence for the occurrence 
of order-disorder pseudo-phase transitions is provided by 
the order parameter and the response functions of finite 
samples. Then, studying the order parameter distribu- 
tion functions, a clearly different behavior between the 
noncritical regular lattice (p = 0) and the small-world 
network [p > 0) is observed. Indeed, the latter shows 
the behavior expected for systems undergoing thermally- 



driven continuous phase transitions. Hence, it is con- 
cluded that a small fraction of shortcuts is sufficient to 
trigger criticality. 

In order to obtain additional evidence of this phe- 
nomenon, standard finite-size scaling relations are used 
to determine the phase diag ram vs p, in which the crit- 
ical temperature is shown as a function of the shortcut- 
adding probability of the network. As expected, for p > 
the system undergoes order-disorder continuous phase 
transitions at finite critical temperatures. Since the long- 
range ordering is weaker the lower the shortcut fraction, 
Tc decreases monotonically with p and vanishes for p — 0. 
Moreover, the behavior of the system at criticality is fur- 
ther characterized by the calculation of the critical expo- 
nents = 3.6 ± 0.4 and 7 = 3.3 ± 0.4. 

These results, obtained in the framework of nonequi- 
librium growth systems, are a novel realization of anal- 
ogous phenomena, which have recently been reported in 
the investigation of related equilibrium systems. Indeed, 
they provide further and complementary evidence on the 
ordering and criticality-inducing effects of long-range in- 
teractions on dynamical systems defined on small-world 
networks. The present work will thus hopefully stimulate 
and contribute to further developments in the fields of 
complex networks and nonequilibrium statistical physics. 
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